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We here report a few numerical tests comparing geometric integrators, of Runge-Kutta type, described by Butcher 
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D. = diag(Z?i, . . . ,bk) and, finally, by considering the Legendre polynomials Pj{i) of degree j — 1, for j > 1, shifted 
and normalized in the interval [0,1] so that /J Pi(x)Pj(x)dz = 5,y (the Kronecker symbol), & = (P/(c,)) g IR SXi , 
3? k — (Pj(cj)) g IR ixi+1 . Method (l)-(a) reduces to the s-stage Gauss-Legendre method when a = (see, e.g., [10, 
pp.77 ff.]). The same happens to method (l)-(b) when k = s [5, 6]. The i-stage Gauss-Legendre method is known 
to be a symplectic integrator of order 2s, able to preserve quadratic invariants for Hamiltonian problems in canonical 
form [9]. On the other hand, under suitable mild assumptions [7] the parameter a in (l)-(a) can be tuned, at each step, 
in order to obtain also the conservation of the Hamiltonian (see also [8]): let us denote such methods by EQUIP(.s') 
(Energy and gf/adratic Invariants Preserving) methods. Finally, the formulae (l)-(b) define the class of HBVM(A:,i) 
methods [2, 3, 4, 5, 6], able to preserve polynomial Hamiltonian functions of degree v, provided that k > (vs)/2 
(obviously, a practical conservation of energy is obtained, for all suitably regular Hamiltonian functions, provided that 
k is large enough). The order of all the above mentioned methods is 2s. In the following we fix s = 3. 

In Figures 1-3 we plot the errors (in the solution, in the Hamiltonian, and in the angular momentum, respectively) 
versus the (constant) stepsize used, for the GAUSS(3) (= HBVM(3,3)), HBVM(4,3), HBVM(6,3), HBVM(9,3), 
HBVM(12,3), and EQUIP(3) methods applied to the Kepler problem [9, pp. 7-9], with eccentricity e = 0.6, over 
1000 periods. As one can see (Figure 1), the order of all methods is confirmed to be 6, even though the error constants 
of HBV(£,3), k > 3, and EQUIP(3) methods turn out to be apparently the same, and approximately 40 times less than 
that of GAUSS(3) (= HBVM(3,3)). The error in the Hamiltonian (Figure 2), as expected, decreases for HBVM(/fc,3) 
methods, as k is increased (with order 2k [4], until round-off errors prevail), and practical conservation is obtained for 
k>9. EQUIP(3) clearly conserves, by its own definition, the Hamiltonian. Finally (Figure 3), the error in the angular 
momentum (which is a quadratic invariant) is negligible for GAUSS(3) and EQUIP(3) methods, and decreasing at the 
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FIGURE 1. 
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FIGURE 2. 



same rate 6 (= 2s) with the stepsize, for HBVM(£,3), k > 3, methods. This is to be expected, since this error only 
depends on matrix A3 (see (l)-(b)), which is the same for all such methods. 

To conclude, we report the numerical results, by using variable stepsize with a standard stepsize selection strategy 
(tol = 1(T 8 ) , for the GAUSS(3), HBVM(12,3), and EQUIP(3) methods applied to the Kepler problem, with eccen- 
tricity e = 0.99, over 100 periods. All methods select stepsizes in the range 10~ 4 ~ 10°. As is well known [9] standard 
stepsize strategies don't work well with symplectic methods, so that GAUSS(3), though preserving the angular mo- 
mentum, exhibits a drift in the numerical Hamiltonian (see Figures 4 and 5). On the contrary, HBVM(12,3) practically 
conserves the Hamiltonian but exhibits a drift in the angular momentum (see Figures 6 and 7). At last, from Figure 8 we 
conclude that only EQUIP(3) preserves both the energy and the angular momentum, when a standard mesh selection 
strategy is used. 
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KEPLER problem, e=0.99, 1 OO periods, tol = 1 O -1 
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KEPLER problem, e=0.99, 1 OO periods, tol = 1 O 
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